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ABSTRACT 

Supersonic  combustion  in  the  inviscid  stream  has  been 
shown  to  reduce  the  base  drag  on  a  projectile.  Fast  reacting 
fuel  burned  externally  in  the  inviscid,  supersonic  flow  gen- 
erates compression  waves,  which,  when  properly  originated, 
impinge  on  the  recirculation  bubble  aft  of  the  projectile 
base  and  increase  the  base  pressure.   A  free- jet  wind  tunnel 
was  constructed  to  study  external  burning  where  heat  addition 
was  simulated  by  nozzle  contour.   Sufficient  increase  in  base 
pressure  was  attained  to  not  only  reduce  base  drag  but  also 
provide  thrust.   Thus,  external  burning  provides  a  method 
for  increasing  the  performance  of  projectiles. 
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I.   INTRODUCTION 

For  an  axisymmetric  object  flying  at  supersonic  speed  as 
much  as  50  percent  of  the  overall  drag  can  be  attributed  to 
the  negative  gage  pressure  acting  upon  the  base.   This  rela- 
tionship between  base  drag  and  the  remaining  drag  on  the  body, 
which  consists  of  wave,  skin  friction,  and  boat-tail  drags, 
is  observed  primarily  in  the  low  supersonic  flight  regime. 
For  Mach  numbers  near  2  the  pressure  on  the  base  can  be  re- 
duced to  as  little  as  50  percent  that  of  the  ambient  air 
while  the  nose  pressure  can  be  kept  relatively  low  by  using 
slender,  streamlined  shapes.   With  increasing  Mach  number, 
however,  a  limit  to  the  base  drag  is  reached  since  the  base 
pressure  cannot  be  less  than  a  negative  one  atmosphere, 
whereas  the  pressure  on  the  nose  will  continue  to  increase 
with  Mach  number. 

The  reduction  of  drag  is  desirable  in  order  to  increase 
the  range  and  terminal  velocity  of  a  projectile.   Pressure 
or  wave  drag  and  skin  friction  drag  have  been  satisfactorily 
reduced  by  designing  streamlined  shapes.   Base  drag,  however, 
has  not  been  reduced  significantly  although  the  use  of  boat- 
tails  and  base  bleed  have  contributed  slightly.   The  most 
promising  method  of  base  drag  alleviation  is  that  due  to 
external  burning.   Heat  addition  in  the  inviscid  flow  in 
proximity  of  the  base  will  create  compression  waves  that, 
if  optimally  originated,  will  impinge  on  the  base  flow 
region.   The  interaction  of  these  compression  waves  will 
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deflect  the  free  streamline  toward  the  symmetric  axis,  and 
the  static  pressure  in  the  recirculation  zone  will  increase 
thus  reducing  the  base  drag. 
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II.   THEORY 


The  base  flow  problem  is  visualized  most  easily  with 
the  inviscid  flow  model  in  Fig.  1  which  ignores  the  effects 
of  viscosity. 

A.   INVISCID  FLOW 

A  semi-infinite,  two-dimensional  model  can  be  formu- 
lated which  easily  satisfies  all  the  necessary  boundary 
conditions. 

1.   Two-Dimensional  Model 

If  a  value  of  P_  (=  base  pressure)  is  selected, 
then  the  initial  turning  through  the  Prandtl  Meyer  expansion 
is  uniquely  determined.   The  pressure  and,  therefore,  the 


EXPANSION    RAN 


8ASE 


PLANE  OF  SYMMETRY 


FIG.    I 

TWO-DIMENSIONAL      INVISCID       FLOW 
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velocity  and  Mach  number  are  constant  along  the  streamline 
BC.   The  Mach  number  can  be  calculated  from  the  Prandtl 
Meyer  equations.   For  uniform,  two-dimensional  flow  over  a 
convex  corner,  the  pressure  depends  only  on  the  angle  of 
inclination  of  the  free  streamline  (angle  0  in  Fig.  1) ,  and 
thus  BC  is  a  straight  line.   The  shock  wave  originating  at 
C  is  weakened  eventually  by  the  interaction  of  the  expansion 
waves.   Near  the  plane  of  symmetry,  the  shock  must  return  the 
flow  to  the  free  stream  direction.   As  the  flow  continues 
downstream  of  the  trailing  shock,  the  static  pressure  ap- 
proaches the  free  stream  pressure.   Unfortunately,  for  a 

given  M   ,  an  infinite  number  of  possible  solutions  exist. 

3      oo  L 

As  PR  is  decreased,  the  Prandtl  Meyer  angle  is  increased, 
and  point  C  moves  closer  to  the  base.   A  limit  is  reached 
eventually  where  for  a  given  local  Mach  number  M. ,  the 
trailing  shock  cannot  return  the  flow  to  its  original 
direction.   The  maximum  deflection  shock  wave  for  a  given 
local  Mach  number  determines  a  limiting  flow  from  the  infi- 
nite number  of  possible  solutions.   This  limiting  case 

yields  the  minimum  possible  PD  and  thus  the  maximum  base 

o 

drag  in  inviscid  flow. 

2.   Axisymmetric  Model 

Inviscid  theory  can  be  extended  to  axisymmetric 
shapes  with  some  complications.   The  expansion  characteris- 
tics issuing  from  the  base  corner  are  not  straight  lines  as 
in  the  Prandtl  Meyer  flow  except  near  the  corner.   Because 
the  flow  conditions  upstream  of  the  trailing  shock  are 
dependent  not  only  on  the  streamline  inclination  at  the 
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base  corner  but  also  on  the  history  of  the  flow  upstream  of 
the  Mach  lines  passing  through  the  corner,  the  free  stream- 
line of  constant  pressure  is  not  straight,  as  shown  in  Fig. 
2.   The  method  of  characteristics  can  be  utilized  in  order 
to  determine  possible  flow  patterns. 


SHOCK 


FIG.  2 
AXI  SYMMETRIC      INVISCID    PLOW 


The  inviscid  flow  pattern  of  an  axisymmetric  body  cannot  be 
constructed  to  the  axis  of  symmetry  and  still  satisfy  the 
boundary  conditions.   This  is  a  consequence  of  the  curvature 
of  the  free  streamline  and  the  fact  that  for  a  given  Mach 
number,  a  single  trailing  shock  is  limited  in  the  amount  of 
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turning  it  can  provide.   The  angle  of  inclination  of  the 
constant  pressure  free  streamline  continues  to  increase 
raonotonically  as  the  axis  of  symmetry  is  approached.   The 
flow  pattern,  therefore,  cannot  be  continued  beyond  a  point 
where  the  limiting  angle  (a  in  Fig.  2)  exceeds  the  angle 
that  the  single  trailing  shock  wave  can  return  the  flow 
parallel  to  the  axis.   In  real  flows,  the  wake  prevents  the 
free  streamline  from  reaching  the  axis,  and  thus  this  limita- 
tion yields  some  similarity  to  observed  flows.   Flow  fields 
containing  a  free  streamline  which  does  not  reach  the  axis 
of  symmetry  can  be  considered  as  those  that  would  exist  in 
an  inviscid  flow  about  an  axisymmetric  body  which  has  an 
infinitely  long  cylindrical  sting  attached  at  the  base  as 
shown  in  Fig.  3. 


M, 


-A2h 


FREE 
STREAMLINE 
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FIG.  3 
CHAPMAN    ANALOGY    BETWEEN    WAKE 
DIAMETER    AND    INVISCID    FLOW    ROD 

DIAMETER 
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The  fraction  d/2h  equals  the  ratio  of  rod  diameter  to 

body  diameter.   As  d/2h  approaches  one,  the  two-dimensional 

result  for  P0  is  obtained.   With  a  rod  present,  inviscid 
is 

theory  predicts  a  lower  base  pressure  and  consequently  a 
higher  base  drag  for  two-dimensional  profiles  than  for 
axisymmetric  shapes. 

This  is  in  qualitative  agreement  with  the  experiment. 
Chapman  [Ref.  1]  found  that  the  base  drag  of  axisymmetric 
bodies  in  an  inviscid  flow  was  so  small  that  it  could  not 
account  for  observed  drag  values  determined  from  experi- 
ment in  real  fluid  flow.   This  strongly  suggested  that 
viscous  effects  were  essential  in  determining  base  flow. 

B.   VISCOUS  FLOW 

The  inclusion  of  viscosity  effects  results  in  a  more 

complicated  flow  model  as  shown  in  Fig.  4.   The  fluid  viscosity 


WAKE 


7p       7*~    yC- 


FIG.  4- 
FLOW    AT    THE.    BASE 
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causes  a  slow  circulatory  motion  in  the  dead  air  region,  and 
the  viscous  mixing  process  causes  the  boundary  layer  to 
thicken  as  it  approaches  the  trailing  shock  wave.   The  base 
pressure  arrives  at  a  steady  state  equilibrium  value  because 
of  the  scavenging  effect  of  the  outside  flow  on  the  recircula- 
tion region.   As  some  of  the  dead  air  mass  is  removed,  the 
angle  of  turning  at  the  base  corner  is  increased,  and  con- 
sequently the  base  pressure  is  decreased.   The  larger  turning 
angle  increases  the  flov;  velocity  outside  the  boundary  layer 
which  increases  the  scavenging  action,  thereby  lowering  the 
base  pressure.   The  process  is  repeated  until  point  C  in 
Fig.  4  reaches  an  equilibrium  position.   As  point  C  moves 
toward  the  base,  however,  the  pressure  ratio  of  the  trailing 
shock  wave  increases,  making  it  more  difficult  for  the 
scavenged  air  and  the  low  velocity  air  in  the  boundary  layer 
to  overcome  the  pressure  rise.  An  equilibrium  position  is 
reached  which  balances  these  two  opposing  effects.  Viscosity, 
therefore,  determines  a  unique  solution  from  the  infinite 
number  possible  for  a  given  Mach  number  and  Reynolds  number 
but,  unfortunately,  complicates  the  theory  excessively. 
C.   CHAPMAN'S  CORRELATION 

Viscous  theory  is  too  involved  to  be  a  practical 
method  for  estimating  base  pressure.   Instead  many  correla- 
tions have  been  attempted  to  aid  in  the  understanding  or  to 
create  a  practical  means  of  predicting  base  drag.   According 
to  Chapman's  model  [Ref.  1]  base  pressure  depends  on  the 
Mach  number  and  pressure  at  the  base  corner,  the  boundary 
layer  thickness  <f ,  the  body  diameter  2h,  and  the  body  angle 
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S  at  the  base.   For  finite  bodies,  as  shown  in  Fig.  5, 
profile  shape  directly  affects  the  Mach  number  and  pressure 
at  the  base  corner  since  the  existence  of  a  bow  shock  wave, 
for  example,  would  alter  the  Mach  number  and  pressure  prior 
to  the  corner. 


BOW 
SHOCK 


EXPANSION 


TRAILIN4G 
SHOCK 


/ 


FIG.   5 
FINITE     BODY     EFFECTS 


Mnand  P.  are  needed  instead  of  M   and  P   .   Chapman 
11  oo      oo 

offered  a  simple  means  of  calculating  the  corrected  approach 
Mach  number  M.  and  pressure  P, .   He  suggested  extending  the 
after  body  one  diameter  (2h)  and  taking  the  average  Mach 
number  and  pressure  along  this  extension  as  the  corrected 
M,  and  P..  .   From  dimensional  analysis,  Chapman  found  that 
C   -  ^C^ij^/h,  jB   )  where  C  R  is  the  base  pressure  coeffi- 
cient and  is  equal  to  (Pg-P-,)  /  ( /?  V,a/2  ).   Viscous  effects 
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were  included  in  a  Reynolds  number  dependency  of  the  boundary 
layer  thickness  <f.   Reynolds  number,  Re  =  /^V^L/Zl^,  ,  is  based 
upon  the  finite  body  length  L  and  the  free  stream  conditions. 
Foo    is  the  free  stream  air  density,  Vo©  is  the  velocity,  and 
Ac?  is  the  coefficient  of  viscosity.   Assuming  a  fixed  M^, 
boundary  layer  theory  predicts  boundary  layer  thickness  £    is 
proportional  to  l/WRe  for  laminar  flow  or  £   is  proportional 
to  l/£Re) s    f°r  turbulent  flow.   With  yS  =  0,  the  final  result 
from  dimensional  analysis  is  C  R  =  J^Ly/^Qe* )  for  laminar  flow 
and  C   =  §(L/\\Rq's   )  for  turbulent  flow. 

Chapman  used  data  from  axisymmetric  flow  to  make  a  corre- 
lation of  C  _  with  L/KRe2,  for  laminar  flow,  as  shown  in 
Fig.  6,  and  L/KRe  s  )    for  turbulent  flow. 

The  effect  of  both  Re  and  body  shape  is  apparent,  and  the 
laminar  correlation  is  reasonably  close  to  the  data.   For 
turbulent  flow,  Chapman's  correlation  again  is  good;  however, 
Re  has  a  much  smaller  effect  on  base  pressure.   This  is  con- 
sistent with  theoretical  considerations  since  l/Re^  ls  much 
less  Re  dependent  than  is  '/Re  x  . 
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D.   EFFECTS  OF  SPIN 

The  effect  of  spin,  or  a  body's  revolution  about  its 
axis  of  symmetry,  on  base  pressure  also  has  been  shown  to  be 
considerable.   Unclassified  material  on  this  subject  is 
lacking,  and  apparently  no  purely  theoretical  approach  has 
ever  been  made  successfully.   It  has  been  concluded  instead 
of  analysis,  that  in  order  to  determine  the  effects  of  spin, 
experimental  measurements  on  specific  object  shapes  need  to 
be  made. 

C.  Wieselsberger  [Ref.  9]  showed  in  19  27  that  the  drag 
coefficient  of  bodies  at  subsonic  speeds  could  be  increased 
by  spin.   Eased  on  drag  data  of  a  projectile  and  cylinder  at 
low  subsonic  speeds  with  circumferential  velocities  up  to 
three  times  the  forward  velocity,  Wieselsberger  concluded  that 
the  increase  of  total  drag  resulted  from  increases  in  skin 
friction  and  base  drag  due  to  spin  induced  alterations  of  the 
boundary  layer  velocity  profile.   R.  Lehnert  [Ref.  9]  showed 
that  spin  influenced  the  aerodynamic  characteristics  of  pro- 
jectiles at  supersonic  speeds,  but  comprehensive  information 
concerning  the  effect  of  spin  was  lacking.   H.  H.  Kurzweg 
[Ref.  7]  showed  a  close  relationship  between  boundary  layer 
and  base  pressure,  as  did  Chapman,  which  prompted  a  study  by 
R.  Lehnert  and  S.  Hastings  [Ref.  9]  to  determine  the  base 
pressure  behavior  as  a  functional  of  spin  rate.   Experiments 
were  conducted  at  the  Naval  Ordnance  Laboratory  in  a  super- 
sonic aeroballistics  tunnel.   Cone  cylinders  of  varying 
lengths  were  spun  at  rates  varying  from  0  to  4  3,000  rpm,  and 
base  pressure  was  measured  through  an  orifice  located  on  the 
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sting  immediately  adjacent  to  the  model  base.  Spark-schlieren 
photographs  were  taken  in  order  to  study  the  effect  of  spin 
rate  upon  the  wake. 

In  Fig.  7  base  pressure  ratio  was  plotted  versus  model 
spin  rate  with  the  dimensionless  model  length  as  a  parameter. 
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The  upper  curve,  which  corresponds  to  the  shortest  model, 
shows  a  considerable  decrease  in  base  pressure  ratio  with  in- 
creasing spin  rate.   The  next  lower  curve  also  shows  a  decrease 

in  P„/P   with  increasing  spin  rate.   The  three  lower  curves, 
B  co 

corresponding  to  the  three  longest  models,  show  relatively 
little  change  in  the  base  pressure  ratio  with  changing  spin 
rates.   As  was  determined  by  both  Kurzweg  and  Chapman,  for  a 
given  Mach  number  and  model  configuration,  the  shape  of  the 
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curve  for  base  pressure  versus  Reynolds  number  is  determined 
by  the  state  of  the  boundary  layer. 

Figure  8  shows  base  pressure  ratio  plotted  against 
Reynolds  number  for  Lehnert  and  Hastings'  data. 
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The  Reynolds  number  was  referred  to  free  stream  conditions 
and  the  wetted  model  length.   The  first  curve  shows  the  zero 
spin  data.   The  large  decrease  in  P_  with  increasing  Re  up  to 
3.6  million  is  accompanied  by  the  existence  of  a  laminar 
boundary  layer  along  the  entire  model  length.   Further  in- 
creases in  Re  resulted  in  a  decrease  in  P_,  and  the  transition 
in  the  wake  mixing  zone  appeared  to  move  upstream.   Base 
pressure  reached  its  minimum  value  when  transition  occurred 
at  the  model  base  and  increased  only  slightly  at  higher  Re 
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when  transition  occurred  upstream  of  the  base.   The  curves 
for  various  spin  rates  show  a  considerable  decrease  in  base 
pressure  with  increased  spin  when  laminar  flow  existed  at 
the  base.   This  decrease  in  base  pressure  due  to  spin  was 
accompanied  by  an  upstream  travel  of  transition  in  the  mixing 
zone  as  was  the  case  for  zero  spin.   In  the  transition  region, 
the  base  pressure  is  unaffected  by  spin  when  the  boundary 
layer  is  turbulent  at  the  base.   Increasing  spin  resulted 
in  a  slight  base  pressure  increase.   The  minima  of  the  base 
pressure  versus  Re  curves  move  toward  lower  Re  with  increas- 
ing spin  rate.   It  can  be  concluded  that  spin  destabilizes 
the  boundary  layer  and  thereby  produces  premature  transition. 
Transition  at  zero  spin  was  observed  to  occur  at  nearly  one 
diameter  upstream  of  the  base  while  at  40,000  rpm,  the  transi- 
tion moved  to  almost  two  diameters  upstream. 

Kurzweg  [Ref.  8]  and  Chapman  [Ref.  1]  observed  that  base 
pressure  and  wake  length  were  correlated.   Decreasing  base 
pressure  corresponds  to  decreasing  wake  length  as  shown  in 
Fig.  9.   This  was  confirmed  by  Lehnert  and  Hastings  [Ref.  9] 
and  suggested  a  method  of  measuring  base  pressure  on  free 
flight  projectiles.   Note  the  similarities  of  Fig.  9  with 
Fig.  7. 
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E.   HEAT  ADDITION 

Base  drag  reduction  by  external  combustion  was  studied 
by  W.  C.  Strahle  [Ref .  13] .   Strahle  observed  that  the  injec- 
tion of  mass  or  heat  into  the  viscous  recirculation  bubble 
resulted  in  low  efficiency  since  the  base  pressure  rise  was 
limited  to  the  free  stream  atmospheric  limit.   Combustion  in 
the  inviscid  supersonic  stream,  however,  imposed  no  such 
limit  and,  therefore,  showed  the  most  promise  for  reducing 
the  base  drag   or  for  providing  a  propulsive  mechanism. 
1.   Control  Volume  Analvsis 

A  control  volume  analysis  [Ref.  6]  as  shown  in 
Fig.  10  gives  some  intuitive  feeling  for  heat  addition  in 
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the  supersonic  stream.   Applying  the  momentum  theorem  for 


the  control  volume  in  the  x  direction,  one  has 

CP6-PJA6 =£e/VdA  -jTVcAM  +  jTVp-P„)<U  ^CP-PjdA      (I) 

where  dA  is  the  projected  area  perpendicular  to  the  free 
stream  direction  x.   For  planar  flow  and  no  heat  addition, 
pressure  P  along  streamline  cd  is  determined  by  the  local 
flow  angle  in  accordance  with  the  Prandtl  Meyer  function. 
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BASE     CONTROL    VOLUME 


Figure  11  shows  that  the  static  pressure  has  a  lower  value 
in  the  fast  turning  case  (short  recirculation  region)  as  com- 
pared to  the  slower  turning  case  (long  recirculation  region) . 

The  /  (P— f^dA  will  be  less  for  a  fast  turning  stream- 
line; and,  consequently ,  PD  also  will  be  lower.   Since  by  the 


con 


tinuity  equation  f  fuAA  -  f  t fuAA  r    only  changes  in  momentum 
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thickness  £**  and  displacement  thickness  ^*  can  change  the 
momentum  flux  through  de  and  be,  it  is  the  gage  pressure 
integrated  over  cd  and  de  that  largely  determines  base 
pressure. 
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Compression  waves  emitted  from  the  heat  addition  zone  as 
shown  in  Fig.  12  increase  the  static  pressure  along  the  edge 
of  the  shear  layer,  cd  in  Fig.  11;  and  therefore,  /  (P-fco'JcA 
is  increased  thereby  increasing  P_.   The  pressure  decreases 
along  cd  as  the  flow  turns  toward  the  axis.   By  adding  heat 
in  the  inviscid  region,  the  shear  layer  is  turned  inward 
with  i\6   or  less  loss  of  static  pressure  along  cd .   PB  is  in- 
creased which  results  in  a  base  drag  reduction. 
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2 .   Strahle  Two-Dimensional  Model 

Strahle  [Ref.  13]  restricted  his  heat  addition 
model  to  two-dimensional,  turbulent,  supersonic  flow,  with 
the  heat  addition  zone  aligned  with  the  inviscid  stream 
beyond  the  base  corner  turn  as  shown  in  Fig.  13. 

The  heat  addition  zone  is  rectangular  in  shape  and  is 
assumed  far  enough  removed  from  the  viscous  region  to  pre- 
clude any  heated  streamlines  being  entrained  in  the  viscous 
flow  region,  which  would  result  in  a  performance  loss. 
Strahle  chose  the  mixing  theory  of  Crocco  and  Lees  [Ref.  2] 
to  describe  the  base  flow  because  of  the  method's  success  in 
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predicting  the  base  pressure.   Strahle  found  that  the 
inclusion  of  external  combustion  did  not  alter  the  base 
flow  theory  of  Crocco  and  Lees  but  instead  only  changed  the 
initial  conditions  for  Crocco  -  Lees  governing  differential 
equation  (Eg.  2.30  in  [Ref.  2]). 
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FIG.    13 
STRAHLE  MODEL  FOR  HEAT  ADDITION 
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POINT 


With  heat   addition,    the  angle   of   inclination    0    in  Fig. 

14    is   no   longer  a   function  of  W      alone  but   is   also  a   function 

of  both  W     and  x. 
e 

The  symbol  W   is  the  ratio  of  the  x  component  of  velocity 
to  the  speed  of  sound  a  .   A  saddle  point  singularity  exists 
through  which  the  trajectory  must  pass  in  order  that  the  far 
wake  limit  can  be  reached.   This  forms  the  uniqueness 
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condition;  since  without  it,  an  infinite  number  of  solutions 
would  exist.   This  saddle  point  is  called  a  critical  point 
and  physically  represents  the  transition  point  from  subsonic 
to  supersonic  flow  in  the  viscous  region.   A  unique  solution 
is  forced  by  the  critical  point  regardless  of  the  external 
combustion  occurring  in  the  inviscid  stream.   The  position  of 
the  critical  point,  however,  and  the  influence  on  the  allowable 
initial  condition  will  be  altered  since  0   is  influenced  by 
combustion. 
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STRAHLE'S   LINEAR    THEORY    CONFIGURATION 


Strahle  analyzed  the  combustion  region  by  the  linear  theory 
of  heat  addition  applied  to  a  supersonic  stream.   The  X  and  Y 
coordinates  defined  in  Fig.  14  are  obtained  by  rotating  the 
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axis  through  an  angle  0(©)and  are  aligned  with  the  rectangular 
heat  source.   In  the  uniform  field  with  Mach  number  equal  to 
M_  ,  the  perturbation  velocity  potential  equation  is 

~Pl  0„  +  0>w  =  5  ■  ( y-ik//"B  asx  (2) 

and  the  pressure  coefficient  is  C  =  ~2.^x:/U&.  Upon  inte- 
grating Eq.  (2)  where  f  is  non  zero  only  for  positive  X  and 
Y,  along  the  Y  =  0  line,  the  perturbation  velocity  potential 

is    0x(X,o)  =  -  1/flSjf  f  CK  (*-*'V0a1<lx.  (3) 

With  0Co}  unknown  and  thus  the  orientation  of  cj_  ,  which  is  the 
heat  release  per  unit  volume  per  unit  time,  also  unknown,  the 
Bernoulli  equation  along  a  streamline  states 

JP/ftf  *-</«/«  -  Ao/tfctt - 60x/u  =  -A(tAy2t*Li*ai-t)nx/il  (4) 

Solving  for<j0  and  integrating  yields 

e-ek=o(Mbvo(MH/?B^/u8  (5) 

where  0  is  the  Prandtl  Meyer  function.  0   is  evaluated  at 
0(o)    on  the  Y=0  line  regardless  of  the  history  of  the  viscous 
region  edge.   For  large  heat  addition,  this  is  a  source  of 
error  since  the  initial  direction  of  the  viscous  region  is 
different  from  its  alignment  at  the  conclusion  of  the  in- 
fluence of  heat  addition.   Integrating  Eq.  (3)  and  substi- 
tuting into  Eq.  (5)  gives 

e-e>b=  OCM^-oc^-^goo  fa) 

where   ^  =  (  j-|)  A^/^Q  a%/i  W 

=    X-Xi-^Y;  if     x,  £  x<xa 

=  Pb  l>  i if    x^^  x<x3 
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with 


Xx=Xi+|B8(Y;+|>) 
X3  =  Xi+Wi.|5BY; 

x4  -  xi+w^B(^-  +  b) 


3 .   Strahle  One-Dimensional  Model 

Because  the  linear  theory  contains  contradictions 
of  the  assumptions  at  even  modest  levels  of  heat  addition, 
Strahle  suggested  a  one-dimensional  approach  for  comparison. 
Figure  1.5  shows  the  basic  scheme  under  a  small  angle  assump- 
tion that   JA^/o'x  =©2.-©,    and  dAc/dx  =  0,  ~  © 

v/here  A^  is  the  area  of  the  heat  addition  zone  and  A   is  the 
n  *= 

area  of  the  adiabatic  free  stream. 


FIG.  15 

strahle's  one-dimensional  heat  addition 
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Assuming  the  flow  on  the  upper  side  of  A,  (Fig.  15) 
contains  waves  of  only  one  family,  then  ©a  =  0^  *  ^t  Mb)- ^(^ )  • 
Using  the  one-dimensional  flow  relations  that  (dAe/Ae^d  %  = 
0-Me)W^P«/Pe)d*'  Strahle  obtains 

©-©b-0(MbV^C^AVCAe/PG)[(l-He1)^Mc]JPe/<jx--d4^/clK.  (7) 

dPe/dx  is  always  maintained  small  by  the  combustion-affected 

region.   Thus,  for  large  heat  release,  dA,/,   dominates  the 
pressure  gradient  term.   If  f^qCx)   is  identified  with  dA,/dx 
in  Eq.  (6) ,  then  an  analogy  exists  between  the  one-dimensional 
and  linear  theory  approaches.   On  the  average,  dA^ /dx  = 
b(AV/v)/W  where  /W/V/  is  the  constant  pressure  volume  expan- 
sion of  the  combustible  mixture,  b  is  the  heat  addition  zone 
capture  height,  and  W  is  the  length  of  the  combustion  zone. 
With  qCx)   assumed  in  the  strongest  region  where  X,  f=  X  ^  Xx 
in  Eq.  (6) ,  the  equivalence  is  made  that 

AV/V  =  (3  ~  0  AKf/A  G»i       a*&  a  =  &V/vWy00  (8) 

where  A  \r\*         is  the  heat  release  per  unit  mass  of  fuel,  AC 
is  the  air  to  fuel  ratio,  and  jS,   s>vMj!|-l  .   Functionally, 
&V/V  -  j(/4.f  Pq,~Tq    )   with  the  most  important  dependences 
being  on  A  and  T~,  which  is  the  base  temperature.   Strahle 
hoped  that  by  combining  linear  theory  with  a  one-dimensional 
approach  the  features  of  both  methods  could  be  included  while 
removing  most  of  the  objections  of  each  method  individually. 
4 .   Strahle *s  Performance  Parameters 

Strahle  computed  the  force  on  the  base  to  be 
Fq   =2(Pb—  F^  )|^   where  P0  is  the  pressure  on  the  base  (where 
X  =  0)  and  h  is  the  base  half  height.   The  specific  impulse 
is  defined  as  X^p s  Fg  Ayw  where  Yn>    is  the  fuel  flow  rate. 
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For  small  values  of  force  level  on  the  base  and  letting  M.-+00 

lip-  PAH, (Jr-0/Me>8B  kWj  (Pe/Pb)fJP/dn,)  (9) 

where  P  =  PB/Pb  and  Kj,,  =  a  heating  parameter.   The  fact  that 
/{•+<*>  creates  an  absolute  linear  heat  addition  limit,  which 
though  impractical,  does  suggest  the  strong  role  that  the 
critical  point  plays  in  the  nonlinear  regime. 

Figure  16  shows  the  effect  the  heat  addition  zone  location 
has  on  specific  impulse.   As  X.  is  moved  from  the  base  toward 
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STRAHLES   SMALL  DISTURBANCE    SPECIFIC  IMPULSE  dP/Ay^ 

VS.  COMBUSTION   ZONE  LOCATION  X; 

the  critical  point,  which  is  located  at  approximately  x  =  7.7 
in  the  linear  limit,  the  specific  impulse  C<JP/dfyJ  reaches  a 
maximum  when  the  combustion  zone  is  roughly  centered  between 
the  base  and  the  critical  point.   A  steep  drop  in  effective- 
ness occurs  as  the  zone  nears  the  critical  point  and  extends 
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beyond  it.   This  is  expected  because  any  influence  upon  the 
reattachment  process  downstream  of  the  critical  point  where 
the  flow  is  supersonic  cannot  be  transmitted  upstream.   The 
linear  limit  predicts  that  external  combustion  will  result 
in  a  base  pressure  rise,  the  specific  impulse  is  expected  to 
be  high,  and  the  combustion  is  wasted  if  the  length  of  the  heat 
addition  zone  W  is  too  long. 

Strahle  recomputed  the  example  in  Fig.  16  for  finite 
combustion  strength,  and  the  resulting  axial  pressure  dis- 
tributions are  shown  in  Fig.  17. 
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FIG.  17 
STRAHLE'S    FINITE  AXIAL    PRESSURE  DISTRIBUTIONS 


Four  of  the  curves  are  for  X.  =  2.0  which  corresponds  to 

the  maximum  I    location  in  Fig.  16.   The  top  curve  is  for 
sp 

Xi  =  0.2.   Apparently,  for  this  case,  it  is  best  to  have  X. 
small.   The  reason  for  this  lies  in  the  nonlinear  behavior 


36 


of  the  critical  point.   In  all  cases,  the  critical  point  be- 
comes entrained  in  the  expansion  portion  of  the  combustion- 
affected  zone.   The  combustion  length  W,  therefore,  fixes  the 
interaction  length,  and  Strahle  suspected  that,  the  most  effi- 
cient placement  of  the  combustion  region  would  be  as  close  to 
the  base  as  possible.   He  also  reasoned  that  the  base  pressure 
increase  resulted  from  the  entrapment  of  the  subsonic  recircula- 
tion bubble  by  compression  waves  emitted  by  the  combustion  zone 
The  singularity  can  occur  only  when  the  inviscid  flow  is 
directed  inward  or  f or 0«£O  .   The  viscous  zone  thickness, 
however,  must  remain  finite.   Since  the  primary  effect  of 
the  combustion  region  is  the  downward  turning  of  the  flow,  the 
last  available  position  for  the  critical  point  appears  to  be 
within  the  combustion-affected  portion  of  the  flow.   The 
flow  passing  the  base  corner  cannot  turn  as  far  without  com- 
bustion, which  in  turn  indicates  a  rise  in  base  pressure. 
5 .   Axisymmetric  Heat  Addition 

Strahle' s  two-dimensional  heat  addition  requires 
modification  for  axisymmetric  application.   A  method  by  Fuhs 
[Ref.  6]  suggests  that  the  change  in  base  drag  resulting  from 
external  burning  is  £*.D  =  &PnAR  where  A   is  the  area  of  the 
base  for  both  planar  and  axisymmetric  bodies  as  shown  in 
Fig.  18. 

Equal  amounts  of  heat  are  released  in  both  the  planar 
and  axisymmetric  cases.   Heat  addition  in  the  planar  case  is 
represented  by  the  lines  above  and  below  the  base  area  and 
in  the  axisymmetric  case  by  a  concentric  circle  of  radius  r_. 
Since  in  both  cases,  the  length  of  heat  addition  zone  is 
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2-TTf"2  t    the  base  area  influenced  by  heat  release  is  in  the 
ratio 

axisymmetric  base   area   _  v-rx/  -r/ir  nri\ 

planar  base   area  'i  /2.TT  r,rx   ~r,/2rx  C«0) 
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FIG.   18 
COMPARISON    OF    BASE  AREA    FOR  PLANAR  AND 

AXISYMMETRIC    GEOMETRY 


The  axisymnietric  geometry  apparently  has  less  base  area 
influenced  by  heat  addition.   Flow  through  an  annular  heat 
addition  region  was  investigated.   It  was  found  that  the 
compression  waves  merged  sooner  with  axial  symmetry,  and  that 
in  contrast  with  the  planar  case,  there  is  an  optimal  r2/r.. 
ratio. 
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III.   APPARATUS 

In  order  to  study  the  effects  of  heat  addition  on  the  base 
pressure  of  an  axisymmetric  body,  a  free  jet  wind  tunnel  was 
constructed  at  the  Naval  Postgraduate  School's  blow  down 
laboratory.   The  test  model  simulated  a  5 -inch  projectile  at 
23,000  feet  and  a  Mach  number  of  2.0,  which  are  typical  condi- 
tions for  a  maximum  range  trajectory.   Similarity  parameters 
were  matched  by  using  a  Mach  2.0  free  jet  coaxial  nozzle 
shown  in  Fig.  19  which  gave  a  Reynolds  number  based  on  radius 
of  1.51  x  106. 


FIG.  19 
COAXIAL     CYLINDER/ NOZZLE 
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A  turbulent  boundary  layer  was  assumed,  and  the  displace- 
ment thickness^*  was  computed  to  be  0.041  inches  at  the 
projectile  base.   This  provides  a  similarity  parameter  of 
displacement  thickness  to  radius  £*/h   of  0.0164. 

The  equipment  consists  basically  of  a  300  psig  air  supply, 
a  glebe  valve,  a  pressure  reducer  valve,  a  90°  elbow,  a  center 
cylinder  support/flow  straightening  section,  a  nozzle  and  an 
overhead  exhaust  hood.   A  4 -inch  pipe  feeds  the  air  from  the 
supply  tank  to  the  on/off  globe  valve.   The  pressure  reducer 
valve  is  followed  by  a  90°  elbow  which  directs  the  air  flow 
vertically  and  provides  a  feedback  pressure  for  the  regula- 
tion of  the  pressure  reducer  valve.   A  2 -inch  spacer  is 
located  between  the  elbow  and  the  center  cylinder  support 
section  and  provides  internal  access  to  the  centerline  cylinder. 
The  centerline  cylinder  is  supported  inside  the  center  cylinder 
support  section  by  flow  straighteners  located  at  the  top  and 
bottom  of  the  3 -foot  section  of  4 -inch  pipe.   A  section  of  6- 
inch  pipe  houses  the  nozzle  on  top  of  the  cylinder  support 
section  and  provides  armor  protection  and  a  support  for  the 
pressure  probe  as  shown  in  Fig.  20.   Figure  21  shows  the 
exhaust  hood  above  the  apparatus,  the  mercury  manometer  board 
used  to  obtain  pressure  data,  the  mercury  vapor  light  source 
and  one  mirror  used  in  the  schlieren  system  for  flow 
visualization. 

Figure  2  2  shows  the  cylinder  extending  from  the  cylinder 
support  section  with  the  nozzle  and  armor  section  removed, 
along  with  the  pressure  gage  connected  to  the  air  supply 
reservoir. 
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FIG.  Zl 
PHOTOGRAPH    OF  TEST  APPARATUS 


42 


FIG.  22 
TEST    CYLINDER 
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The  cylinder  base  is  formed  by  a  cap  which  contains  six 
pressure  taps.   Four  taps  are  located  across  the  base  radius, 
and  two  taps  are  located  just  upstream  of  the  base  corner. 

The  pressure  from  each  tap  is  carried  downward  through 
the  cylinder  and  the  cylinder  support  in  plastic  tubes  which 
exit  the  apparatus  through  the  2 -inch  spacer  located  on  top  of 
the  elbow.   The  tubes  are  connected  to  a  mercury  manometer  from 
which  the  pressure  data  was  read.   Total  and  static  pressure 
also  can  be  measured  by  the  external  pressure  probe,  which 
will  be  the  primary  base  pressure  source  when  cylinder  spin 
is  included.   The  results  of  this  study  were  for  the  static 
case  with  all  similarity  parameters  matched  except  for  spin. 
However,  because  spin  appears  to  have  a  considerable  impor- 
tance on  base  pressure,  provisions  were  included  to  allow 
spinning  the  projectile.   V7hen  a  5-inch  projectile  is  fired, 
the  barrel  rifling  spins  the  projectile  to  250  revolutions 
per  second.   This  establishes  a  spin  similarity  parameter  of 
r^/u  =  0.1595  where  r  is  the  radius  of  the  projectile,^  is 
the  angular  velocity  and  u.  is  the  free  stream  velocity  at 
Mach  2.0  at  23,000  feet  altitude.   The  similarity  parameter 
will  be  matched  in  the  free  jet  experiment  by  spinning  the 
cylinder  at  4  8,000  rpm.   The  spinning  cylinder  design  con- 
sists of  a  6.5-inch  stainless  steel  cylinder  which  slides 
over  a  support  cylinder  containing  high  pressure  nitrogen 
orifices  to  provide  gas  bearing  support  and  align  the 
spinning  cylinder  radially.   A  small  internal  impulse  turbine 
driven  by  high  pressure  nitrogen  will  spin  the  cylinder. 
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Pleat  addition  was  simulated  by  nozzle  contour.   As  de- 
scribed in  Ref.  6,  the  external  burning  process  can  be 
simulated  by  nozzle  wall  contours  which  generate  character- 
istics identical  to  the  supersonic  combustion  of  an  external 
burning  assisted  projectile.   A  gas  generator  within  the 
projectile  injects  fuel-rich,  hot,  combustible  mixtures  into 
the  supersonic  inviscid  flow  where  ignition  occurs  due  to 
high  temperature  and  the  chemical  nature  of  the  fuel.   The 
region  of  heat  addition  was  assumed  axisymmetric  about  the 
cylindrical  body  of  the  projectile  with  the  radius  r_  in 
Fig.  18  equal  to  the  body  diameter.   For  stoichiometric  fuel 
mixtures,  calculations  showed  that  the  combustion  would  result 
in  a  4.5°  deflection  of  the  streamline  bounding  the  heat 
addition  zone,  and  it  was  this  4.5°  deflection  inward  of  the 
streamline  at  r~  from  the  cylinder  centerline  which  governed 
the  wall  contour.   Half  stoichiometric  nozzles  also  are 
planned  for  future  analysis.   An  iterative  procedure  using 
the  method  of  characteristics  was  used  to  design  the  nozzles 
to  give  the  proper  -4.5°  flow  angle.   Using  W  to  denote  the 
length  of  the  combustion  zone  (see  Fig.  14)  and  h,  the  cylinder 
radius,  several  nozzles  of  varying  V7/h  were  fabricated  with  an 
epoxy  cast  poured  over  a  wax  plug  to  the   correct  contour. 
These  nozzles  provide  several  data  points  shown  in  Table  I 
which  show  the  effect  external  heating  has  on  base  pressure. 
During  testing,  spacers  were  used  to  vary  the  axial  location 
of  the  heat  addition  zone  with  respect  to  the  base  corner. 
Figure  23  is  a  photograph  of  the  nozzles,  spacers,  and  tunnel 
armor  section. 
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TABLE    1 
TEST    MATRIX 
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The  symbol  1  represents  the  distance  from  the  corner 
where  the  first  fully  developed  compression  wave  impinges 
the  cylinder  wall  or  the  wake  downstream  of  the  base,  and  h 
again  is  the  cylinder  radius  as  shown  at  the  top  of  Figs. 
40,  41,  and  42.   Five  1/h  values  were  observed  for  each 
nozzle  as  shown  in  Table  I.   Negative  1/h  denotes  compression 
beginning  upstream  of  the  base  corner,  positive  1/h  signifies 
compression  beginning  downstream  of  the  corner,  and  0  denotes 
compression  starting  at  the  corner.   Table  I  is  the  test 
matrix  obtained  to  date  showing  the  results  of  the  completed 
tests. 
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IV.   TEST  PROCEDURE 

After  the  proper  nozzle  and  spacer  combination  was 
installed  on  the  apparatus  and  the  air  reservoir  was  charged 
to  3  00  psig,  a  data  point  for  the  test  matrix  was  taken.  The 
procedure  involved  opening  the  globe  valve  and  setting  the 
proper  output  pressure  from  the  pressure  reducer  valve.   A 
pressure  tap  in  the  elbow  was  connected  to  a  separate 
pressure  gage  which  was  used  to  set  the  conditions  for  each 
run.   This  tap  also  provided  feedback  pressure  for  the 
pressure  reducer  valve.   With  the  proper  output  pressure  set, 
a  schl.ieren  concave  mirror  system  was  used  to  photograph 
each  test  run  as  shown  in  Figs.  25-39.   Polaroid  positive/ 
negative  film  was  used  in  order  to  facilitate  print 
reproduction.   Base  pressure  information  was  read  from  a 
mercury  manometer  and  recorded  for  data  reduction.   At  the 
conclusion  of  the  data  readings,  the  pressure  reducer  valve 
was  turned  off,  and  the  air  reservoir  was  recharged  for  the 
next  run. 
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V.   RESULTS  AND  CONCLUSIONS 
The  pressure  ratios  (PD/P_ )  as  shown  in  Table  I  indicate 

D  GO 

the  performance  of  the  observed  simulated  external  burning. 
Figure  24  is  a  plot  of  the  results. 
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The  data  points  taken  for  the  two  uniform  flow  nozzles 

compare  favorably  with  M  =  2  results  compiled  by  Przirembel 

and  Page  [Ref.  10],   These  points  represent  the  base  pressure 

ratio  of  the  cylinder  with  no  heat  addition;  and  because 

P-./P   is  less  than  1.0,  base  drag  results.   All  other  data 
B   oo 

points,  except  one,  show  an  increase  in  base  pressure  ratio 
with  the  applied  combustion,  indicating  a  reduction  in  base 
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drag.   All  the  data  for  W/h  =  4.34,  and  the  last  two  points 

for  W/h  =  1.84  in  fact,  indicate  thrust  developed  at  the  base 

since  P^/P    for  these  cases  is  greater  than  1.0.   The  data 

indicate,  as  expected,  that  a  longer  heat  addition  zone 

(Larger  W)  will  result  in  a  higher  P„/P    increase,  since  a 

r      B   oo 

longer  interaction  length  of  compression  waves  acting  upon 
the  recirculation  bubble  has  more  effect  than  a  shorter  zone. 
It  also  appears  that  maximum  benefit  from  the  heat  addition 
occurs  when  the  entire  recirculation  bubble  is  embedded  in 
the  combustion-created  compression  waves.   For  W/h  =  4.34, 
a  decline  in  performance  is  noted  as  1/h  exceeds  -1.0,  which 
indicates  that  for  a  long  heat  addition  zone,  optimal  axial 
placement  is  necessary  so  as  to  not  waste  compression  waves 
downstream  of  the  sonic  point.   Because  supersonic  flow  exists 
downstream  of  the  sonic  point,  any  compression  effect  occurring 
beyond  that  point  will  not  be  transmitted  upstream  into  the 
subsonic  circulation  bubble,  which  would  be  necessary  for  a 
base  pressure  increase.   For  the  shorter  heat  addition  zones 
(W/h  =  1.84  and  0.72),  the  base  pressure  ratio  increases  as 
the  zone  is  moved  rearward.   It  is  assumed  that  no  compression 
was  wasted  beyond  the  sonic  point  up  to  1/h  =  0.5.   Although 
not  tested,  a  decline  in  performance  is  expected  with  further 
rearward  movement  of  the  heat  addition  region  with  respect  to 
the  cylinder  base  corner.   Figures  25-39  show  a  schlieren- 
derived  photograph  of  each  data  point  obtained  with  a 
pictorial  description  of  the  heat  distribution  actually 
simulated.   The  double  cross-hatched  regions  depict  anomalous 
compression,  and  the  single  cross-hatched  areas  show  the  heat 

addition  zones  without  multiple  reflection. 
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FIG.  25 
PHOTOGRAPH       AMD    SCHEMATIC     FOR  w/h-  434,  Jf/M-2.23 
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FIG.  26 
PHOTOGRAPH    AND   SCHEMATIC     FOR  W/h=  4.34,  .?/k=- 0.95 
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FIG.  28 
PHOTOGRAPH   AMD    SCHEMATIC    FOB   w/^ 4.34,^/^  =  0.0 
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FIG.  29 
PHOTOGRAPH  AND    SCHEMATIC    FOR    W/K=  4.34, J?/M  O.SO 
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PHOTOGRAPH  AND      SCHEMATIC     OF   W/h=  1.84,  Jt/hs -2.23 
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PHOTOGRAPH    AND    SCHEMATIC     FOR  W/h  =  |  ©4,jP/h  =  -0.95 
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PHOTOGRAPH    AND    SCHEMATIC     FOR  W/K  =  1. 84,  Jl/k-O.so 
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FIG.  33 
PHOTOGRAPH  AND    SCHEMATIC     FOR  W/^=  I.  84-,  J?/h  =  O.O 


60 


=  1.074 


OO 


Fl  G.  34- 
PHOTOGRAPH   AND   SCHEMATIC     FOR   W/K=  1.84,  Sl/Vs  -  O.50 
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FIG.     36 
PHOTOGRAPH  AND    SCHEMATIC    FOR   W/h  -0.7Z,  J?/h  =  -0.9S 
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FIG.  37 
PHOTOGRAPH    AND    SCHEMATIC     FO&    \AJ/^0.72,$/h*  -O.SQ 
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FIG.  38 
PHOTOGRAPH    AND  SCHEMATIC    FOR  W/h=0.72,  JP/K=  O.O 
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FIG.  39 
PHOTOGRAPH  AND     SCHEMATIC     FOR   VJ/\r\  -0.7Z,  J?A  =  0.  SO 
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A  study  of  the  pictures  shows  that  a  uniformly  distrib- 
uted heat  source  was  not  obtained  as  desired.   Instead,  the 
starting  compression  waves,  which  originated  immediately 
downstream  of  the  nozzle  throat,  reflected  from  the  cylinder 
wall  to  the  nozzle  corner  and  then  into  the  base  flow  region 
as  shownat  the  top  of  Figs.  40,  41,  and  42.   For  the  W/h  = 
4.34  nozzle,  the  simulated  heating  distribution  is  shown  in 
the  middle  of  Fig.  40. 

Compression,  originating  from  the  simulated  heat  addition 
zone,  starts  at  point  b  and  increases  to  uniform  strength  at 
point  e  where  it  levels  at  a  constant  value  until  just  prior 
to  point  h.   The  reflected  compression  waves  from  the  nozzle 
corner  superimposed  upon  the  uniform  heat  cause  a  spike  which 
is  approximately  twice  the  magnitude  of  the  intended  uniform 
value.   This  spike  is  seen  in  the  W/h  =  4.34  photographs  as 
a  shock  wave,  formed  by  a  coalescing  of  the  reflected 
compression  waves,  as  seen  in  Figs.  25-29. 

Figures  41  and  4  2  show  the  expected  compression  wave  re- 
flection pattern  and  heat  distribution  for  W/h  =  1.84  and 
W/h  =  0.7  2  respectively.   Basically,  the  same  phenomena  of 
compression  wave  reflection  occurred.   The  shorter  heat  addi- 
tion zones,  however,  ended  prior  to  the  reflection,  and, 
therefore,  the  reflected  waves  created  an  almost  equal 
strength  combustion  region  downstream  of  the  intended  heat 
addition  zone  rather  than  a  spike  as  seen  for  the  W/h  =  4.34 
case. 
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The  photographs  show  that  for  the  W/h  =  1.8  4  case,  the 
reflected  compression  waves  never  coalesce  into  a  shock  but 
are  spread  over  a  slightly  longer  axial  distance.   For  W/h  = 
0.72,  a  shock  did  form,  but  its  strength  is  notably  less  than 
the  shock  occurring  for  W/h  =  4.34. 

Although  the  reflected  compression  waves  existed  for  each 
data  point,  it  is  assumed  that  the  base  pressure  was  affected 
by  the  reflection  for  only  the  points  circled  in  Fig.  24. 
The  circled  points  were  determined  by  the  location  where  the 
reflected  wave  impinged  the  recirculation  bubble.   As  pre- 
viously described,  base  pressure  can  be  influenced  by 
compression  waves  only  if  the  waves  interact  with  the  base 
flow  prior  to  the  critical  point.   Any  interaction  downstream 
of  the  critical  point  will  not  influence  pressure  at  the  base, 
The  bottom  of  Figs.  40,  41,  and  4  2  shows  the  estimated 
critical  point  location  for  the  three  W/h  nozzles. 

A  vertical  line  extended  through  the  intersection  of  a 
given  1/h  value  and  the  critical  point  curve  will  have  on 
its  left  the  portion  of  the  heat  addition  region  affecting 
the  base  pressure. 

In  conclusion,  it  is  evident  that  external  burning  is 
capable  of  not  only  reducing  base  drag  but  also  producing 
thrust  under  the  proper  conditions.   The  optimal  benefit  is 
derived  when  the  entire  recirculation  bubble  is  embedded  in 
combustion  created  compression  waves.   Waves  beginning  up- 
stream of  the  base  corner  have  limited  value  for  raising  base 
pressure,  and  waves  interacting  downstream  of  the  critical 
point  are  wasted. 

71 


LIST  OF  REFERENCES 


1.  Chapman,  D.  R. ,  "An  Analysis  of  Base  Pressure  at  Super- 
sonic Velocities  and  Comparison  with  Experiment",  NACA  TN 
2137,  July  1950. 

2.  Crocco,  L.  and  Lees,  L.  ,  "A  Mixing  Theory  for  the  Inter- 
action Between  Dissipative  Flows  and  Nearly  Isentropic 
Streams",  Journal  of  Aeronautical  Sciences,  Vol.  19,  No. 
10,  October  1952. 

3.  Darling,  J.  A.  and  Knott.,  J.,  "Effect  of  Spin  on  the  Aero- 
dynamics of  the  20mm  Projectile",  Naval  Ordnance  Laboratory, 
TR  71-180,  September  1971. 

4.  Davis,  Dale  D. ,  "Extension  of  Simplified  Mixing  Theory  to 
Axially  Symmetric  Supersonic  Uake  Flows  and  Application  to 
the  Base  Pressure  Problem  for  a  Body  of  Revolution", 
Princeton  University  Master  of  Science  Thesis  24-V-1952. 

5.  Dickinson,  Elizabeth  R. ,  "Some  Aerodynamic  Effects  of 
Varying  Body  Length  of  a  Spinning  Projectile",  Ballistic 
Research  Laboratories,  MR  1664,  July  1965. 

6.  Fuhs,  A.  E.,  "Quasi  Area  Rule  for  Heat  Addition  in  Transonic 
and  Supersonic  Flight  Regime",  AFAPL-TR-7  2-10 ,  /August  1972. 

7.  Fuhs,  A.  E.,  Smithey,  W.  ,  and  Naber,  M.  ,  "Axisymmetric 
Analysis  of  External  Burning  Assisted  Projectile",  9th 
JANNAF  Combustion  Meeting,  CPIA  Publication  231,  Vol.  II, 
December  1972,  pp.  389-416. 

8.  Kurzweg,  H.  H.,  "Interrelationship  Between  Boundary  Layer 
and  Base  Pressure",  Journal  of  the  Aeronautical  Sciences, 
Vol.  18,  No.  11,  November  1931. 

9.  Lehnert  R. ,  and  Hastings,  S.  M. ,  "Spin  Effects  on  Base 
Pressure  of  Cone  Cylinders  at  M  =  2.86",  Naval  Ordnance 
Laboratory,  TR  2956,  October  1953. 

10.  Przirembel,  G.  E.,  and  Page,  R.  H.,  "Analysis  of  Axisym- 
metric Supersonic  Turbulent  Base  Flow" ,  Heat  Transfer  and 
Fluid  Mechanics  Institute  Proceedings,  pp.  258-27  2,  1968. 

11.  Schmidt,  L.  E.,  and  Murphy,  C.  H.,  "Effect  of  Spin  on 
Aerodynamic  Properties  of  Bodies  of  Revolution",  Ballistic 
Research  Laboratories  Report  715,  November  1953. 

12.  Sedney,  R. ,  "Review  of  Base  Drag",  U.  S.  Army  Ballistic 
Research  Laboratories  Report  1337,  October  1966. 


72 


13.   Strahle,  W.  C,  "Theoretical  Consideration  of  Combustion 
Effects  on  Base  Pressure  in  Supersonic  Flight",  Twelfth 
Symposium  on  Combustion,  Combustion  Institute,  Pittsburgh, 
Pennsylvania,  1969. 


73 


INITIAL  DISTRIBUTION  LIST 


No.  of  Copies 


1.  Defense  Documentation  Center  2 
Cameron  Station 

Alexandria,  Virginia   22314 

2.  Library,  Code  0212  2 
Naval  Postgraduate  School 

Monterey,  California   93940 

3.  Professor  A.  E.  Fuhs  5 
Department  of  Aeronautics 

Naval  Postgraduate  School 
Monterey,  California   93940 

4.  Ltjg.  Gary  J.  Caswell,  USN  3 
399  Troy 

Aurora,  Colorado   80010 

5.  Dr.  Warren  Strahle  1 
Aerospace  Engineering 

Georgia  Institute  of  Technology 
Atlanta,  Georgia   30332 

6.  Mr.  Dominic  J.  Monetta  2 
Gun  Systems  Engineering 

Naval  Ordnance  Station 
Indian  Head,  Maryland   206  4  0 

7.  Dr.  Alan  Roberts  2 
Gun  Systems  Engineering 

Naval  Ordnance  Station 
Indian  Head,  Maryland   20640 

8.  Dr.  Frederick  Billig  1 
Applied  Physics  Laboratory 

Johns  Hopkins  University 

8621  Georgia  Avenue 

Silver  Springs,  Maryland   20910 

9.  Professor  Bruce  Reese  1 
Chairman  of  Aeronautical  Engineering 

Purdue  University 
Lafayette,  Indiana   47907 

10.   Cmdr.  William  Smithey  3 

SMC,  Box  23  4  8 
Naval  Postgraduate  School 
Monterey,  California   93940 


74 


11.  Dr.  James  S.  Holdhusen 
Vice  President 

FluiDyne  Engineering  Corporation 
5900  01 sen  Highway 
Minneapolis,  Minnesota   554  22 

12.  Mr.  Thomas  Curran 

Chief  of  Rani  jet  Technology  Branch 

Air  Force  Applied  Physics  Laboratory 

Wright  Patterson  Air  Force  Base,  Ohio   45433 

13.  Lt.  Noel  P.  Horn,  USN 
SMC,  Box  1183 

Naval  Postgraduate  School 
Monterey,  California   93940 

14.  Dr.  Grant  Hosack 
Rocketdyne 

Rockwell  International 

6  63  3  Canoga  Avenue 

Canoga  Park,  California   91304 

15.  Mr.  C.  Miner 

U.  S.  Army  Missle  Command 
Redstone  Arsenal,  Alabama 

16.  Dr.  Richard  Weiss 

Air  Force  Rocket  Propulsion  Laboratory 
Edwards  Air  Force  Base,  California 

17.  Professor  John  Clarke 
Cranfield  Institute  of  Technology 
Bedford,  England 

18.  Dr.  E.  G.  Broadbent 

Royal  Aeronautical  Establishment 
Farnborough ,  England 

19.  Chairman,  Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California   93940 

20.  Dr„  Ing.  G.  Winterfeld 
D  F  V  L  R 

Institut  f.  Luftstrahlantriebe 

505  Porz-Wahn,  Linder  Hone,  Germany 


75 


SECURITY   CLASSIFICATION   Of-    THIS   PAGE   (Wh  :?rr,d> 


REPORT  DOCUMENTATION  PAGE 

R  E  A  D  I N  S  T R  U  C  TIO N  S 
BEFORE  COMPLETING  FORM 

1.     REPORT   NUMBER 

2.  GOVT    ACCESSION  NO. 

3.     RECIPIENT'S  CATALOG  NUMBER 

4.     TITLE  (end  Subtitle) 

Base   Drag   Reduction  by   External   Burning 

5.     TYPE  OF   REPORT   ft   PERIOD  COVERED 

Master's   Thesis 
September   1973 

6.     PERFORMING  ORG.   REPORT  NUMBER 

7.     AUTKOHCs; 

Gary  Joe   Caswell 

6.     CONTRACT  OR  GRANT  NUMBER,'*; 

9.     PERFORMING  ORGANIZATION    NAME   AND   ADDRESS 

Naval   Postgraduate  School 
Monterey,    California      93940 

10.     PROGRAM   ELEMENT,  PROJECT.   TASK 
AREA  ft   WORK  UNIT  NUMBERS 

II.     CONTROLLING  OFFICE   NAME   AND   ADDRESS 

Naval   Postgraduate   School 
Monterey,    California      93940 

12.     REPORT   DATE 

September   197  3 

13.     NUMBER  OF    PAGES 

77 

14.     MONITORING  AGENCY  NAME  ft    ADDRESS!-//  different  from  Controlling  Office) 

Naval    Postgraduate   School 
Monterey,    California      93940 

1S.     SECURITY  CLASS,  (of  thle  report) 

Unclassified 

IS*.     DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

16.     DISTRIBUTION   STATEMENT  (of  this  Report) 

Approved   for   public   release;    distribution  unlimited. 

17.     DISTRIBUTION  STATEMENT  (of  the  abatrect  entered  In  Block  20,  If  different  from  Report) 

18.     SUPPLEMENTARY  NOTES 

19.     KEY  WORDS  (Continue  on  reverse  aide  If  neceaaasy  and  Identify  by  block  number) 

Base   Drag   Reduction 
Axisymmetric  Heat  Addition 
Spin  Effects   on  Base  Pressure 

20.     ABSTRACT  (Continue  on  reverae  tide  If  noceeaary  and  Identify  by  block  numbor) 

Supersonic  combustion   in  the   inviscid    stream  has   been   shown 
to   reduce   the  base   drag  on   a  projectile.      Fast   reacting   fuel 
burned  -externally   in   the    inviscid,    supersonic    flow   generates 
compression  waves,    which,    when   properly   originated,    impinge   on 
the   recirculation   bubble   aft   of   the    projectile   base   and    increase 
the  base   pressure.      A   free- jet  wind   tunnel  was   constructed   to 
study   external   burning  where   heat   addition  was    simulated  by 

DD     FORM 

u«-'     1   JAN   73 

(Page    1) 


1473 


EDITION   OF    I   NOV  65  IS  OBSOLETE 
S/N    0102-014-  6601   | 


76 


SECURITY  CLASSIFICATION   OF   THIS  PAGE  f»7i«n  Dmtm  r.ntared) 


CliCUHlTY   CLASSIFICATION   OF   THIS  P  *.GE(Whor}  Data  Entered) 


20.   ABSTRACT  (Continued) 

nozzle  contour.   Sufficient  increase  in  base  pressure 
was  attained  to  not  only  reduce  base  drag  but  also 
provide  thrust.   Thus,  external  burning  provides  a 
method  for  increasing  the  performance  of  projectiles. 


DD  Form   1473     (BACK) 

1  Jan  73    i  , 

S/N      0102-014-G601  SECURITY   CLASSIFICATION   OF   THIS  PAGE.fKfi&n  Dflfo  Entered) 

77 


Thesis 

C3255     Caswel 

c.l 


146323 


Base  drag  reduction 
by  external    burning. 


Thesis 

C3255 

c.l 


146323 


Caswel 1 

Base  drag   reduction 
by  external    burning. 


thesC3255 

Base  drag  reduction  by  external  burning 
MB"-' 


3  2768  002  09106  8 

DUDLEY  KNOX  LIBRARY 


